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Abstract 

A finite difference formulation is presented 
for sound propagation in a rectangular two-dimen- 
sional duct without steady flow. Before the dif- 
ference equations are formulated, the governing 
Helmholtz equation is first transformed to a form 
whose solution tends not to oscillate along the 
length of the duct. This transformation reduces 
the required number of grid points by an order of 
magnitude. Example solutions indicate that stepped 
noise source profiles have much higher attenuation 
than plane waves in a uniform impedance liner. 

Also, multiple stepped impedance liners are shown 
to have higher attenuation than uniform ducts if 
the impedances are chosen properly. For optimum 
noise reduction with axial variations in impedance, 
the numerical analysis indicates that for a plane 
wave input the resistance should be near zero at 
the entrance of a suppressor duct, while the react- 
ance should be near the optimum value associated 
with the least-attenuated mode in a uniform duct. 

A close-packed array of quarter wave tubes tuned 
to the proper impedance may approximate this opti- 
mum condition in the entrance region of the liner. 
The low resistance entrance region of a segmented 
liner redistributes the acoustic energy nearer to 
the wall where it can be more effectively absorbed. 

Introduction 

In engine ducts which have variations in axial 
impedance, area, or flow, the finite difference 
approach is an attractive alternate to the conven- 
tional analytical techniques currently used in 
noise propagation analysis. In Refs. 1 to 3, the 
conventional finite difference theory has been 
applied to the problem of sound propagation in 
two-dimensional straight hard- and soft-wall ducts 
with and without steady flow. In Ref. 4, the dif- 
ference technique was applied to the variable area 
case. In all the above, the continuous acoustic 
field is lumped into a series of grid points at 
which the acoustic pressure and velocity are de- 
termined (as in Fig. 1). For long ducts and high 
frequency sound propagation, large computer stor- 
age is required. To eliminate this storage re- 
quirement, a wave envelope concept was developed^) 
which greatly reduced the number of grid points 
necessary to solve a problem. However, under cer- 
tain conditions, the author indicated in Ref. 5 
that the transformation was not appropriate or 
advantageous. This paper presents a generalization 
which removes the limitations of the earlier 
theory. 

Herein, a finite difference technique is for- 
mulated for sound propagation in a rectangular 
two-dimensional duct without steady flow. First, 
the governing Helmholtz equation is transformed 
into a form which reduces the required number of 
grid points in the conventional finite difference 
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theory by an order of magnitude. This technique 
is labeled the Generalized Wave Envelope Concept. 
Next, this concept is applied to a duct with vari- 
able axial impedance for a plane wave input. In 
addition, the effect of stepped noise source pro- 
files on the attenuation of uniform impedance 
liners is investigated. 

List of Symbols 

A coefficient matrix 

c velocity of sound 

AdB decrease in decibels 

E acoustic power 

F column vector 

f frequency 

H duct height 

I dimensionless intensity 

i 

L dimensionless length of duct 

to total number of grid rows (points in y-direc- 
tion) 

n total number of grid columns (points in x-di- 

rection) 

P dimensionless pressure, P(x,y), Eq. (1) 

p dimensionless pressure, p(x,y), Eq. (7) 

'p dimensionless column pressure vector 

S source strength 

x dimensionless axial coordinate 

Lx axial grid spacing 

y dimensionless transverse coordinate 

Ay transverse grid spacing 

Z impedance 

£ dimensionless specific acoustic impedance 

H dimensionless frequency, Eq. (2) 

Tl + wave envelope frequency 

X pressure wavelength 

V*" effective wavelength 

Subscripts 

i,j axial and transverse indexes, see Fig. 1 
Superscripts 

(1) real part 

(2) imaginary part 

Governing Equations 

The governing equation describing the steady 
state propagation of sound is the classic Helmholtz 
equation (dimensionless form) 




An example will be presented shortly to illustrate 
how Eq. (7) can remove the oscillatory nature of 
the pressure. 


where the dimensionless frequency n is given as 



( 2 ) 


These and all other symbols used in the report are 
defined in the List of Symbols. 


For plane waves propagating in hard wall ducts, 
is known to be c/f- For this case Eq. (7) 
reduces to 

JL-R + iLX _ i47TT) = o (8) 

3x Z Zy l dx 


For a hard wall duct infinite in extent, the 
one-dimensional solution of Eq. (1) for pressure 
yields (i) 

p p(i) + iP^ 2 ) at cos(2tiT|x) - i sin(2flrix) (3) 

The analytical and finite difference solutions 
for the acoustic pressure of an example case are 
presented in Fig. 2. To accurately describe the 
spatial pressure profiles in Fig. 2, the required 
number of grid points in the x-direction is 

n « 12 n(j~) (4) 

For the example shown in Fig. 2, the dimensionless 
frequency n equals 1 and (L/H) is 1; thus, 12 
points are necessary. 

Next, consider the example case of a two- 
dimensional soft-wall duct with r\ - 1 and 
L/H =3. A typical pressure profile in the 

suppressor duct is shown in Fig. 3 by the heavy 
solid line, while the dashed line represents the 
envelope of the pressure wave amplitude. Ftohi 
Eq* (4), the number of grid points n should be 
36. However, if the Helmholtz Eq. (1) could be 
transformed so that it would describe the envelope 
(dashed line) of the pressure, the grid points 
could be reduced to 3 or less, for this particular 
case as shown in Fig. 3. Since the total number 
of grid points used is the product of the points in 
the y-direction with n points in the x-direction, 
as shown in Fig. 1, this represents a tremendous 
savings in computer storage and operating time. 

The task of transforming the Helmholtz equation to 
the wave envelope equation is now considered. 


General ized Wave Envelope Concept 


To remove some (if not all) of the axially 
oscillatory part of the wave pressure profile, 
assume 


P(x,y) - p (x , y ) e -i2Trri + x ( 5 ) 


where P represents the pressure of the oscilla- 
ting curve in Fig. 3 and p represents the pres- 
sure of the wave envelope shown by the dotted line 
in Fig. 3, and where 



( 6 ) 


with A + representing the effective axial wave- 
length of the pressure in the duct (Fig, 3). Sub- 
stituting Eq. (5) into the Helmholtz equation (1) 
yields a new governing differenital equation called 
the Generalized Wave Envelope Equation 


^ 2 +^ - iW- + (2ir)2( n 2 _ n + 2)p « 0 (7) 
9x Z 


In soft wall ducts, is not known precisely; 
therefore, the problem of picking A + (or n + ) to 
remove the pressure oscillations must be considered. 
In Ref* 5, ^soft was assumed equal to ^hard* 

Later in this section, we will show how the assump- 
tion works. 

For convenience, the hard wall case, for which 
the correct answer is known, will now be used to 
investigate the sensitivity of an incorrect guess on 
the final answer. To illustrate how Eq. (7) or (8) 
removes the oscillatory nature of the pressure, 
consider again the hard wall duct problem of Fig, 2. 
The exact analytical solution of Eq. (7) for p 
in which the assumed n + term is included is 

p or p(l) + ip (2) * cos27r(n-n + )x - isin27r(n-r) + ) x (9) 

Of course, the product of Eq. (9) with e~i 2lT n +x 
always gives the exact pressure variation P as 
defined by Eq. (3). That is, every p(x) obtained 
for any assumed value of n+ will yield the cor- 
rect value of F(x) when substituted back into 
Eq. (3). 

The analytical solutions for pO.) from Eq. (9) 
for various assumed values of n + or |n”Tl + | are 
shown in Fig. 4. For n + equals 0, Eq. (7) reduces 
to the Helmholtz equation, and Pw = p(l), see 
dotted line (cosine shape) in Fig. 4. For this 
particular problem twelve grid points would be re- 
quired to adequately describe this pressure profile 
in a difference analysis. However, if n + is as- 
sumed to have a value of q + between 0.7 and 1.3 
such that |n“n + |<0.3, the exact analytical solution 
for p(l) would require as few as 3 grid points in 
the difference analysis. The curves for n + be- 
tween 1,0 and 1.3 are mirrow image to those between 
0.7 and 1.0 and therefore are not shown. If is 
assumed to be equal to f|, the curve for p(I) is 
the straight line, n = 1 or jrm + [ = 0. This 
line represents the envelope of the pressure oscil- 
lation. If, for example, it is assumed that 
ri + 65 .8, the exact solution in Fig. 4 does not rep- 
resent the true envelope of the pressure oscilla- 
tion; it is still a gently varying function that 
requires only a few finite difference grid points 
to accurately describe its shape. Thus, it is only 
necessary to pick a value of X+ (or rf*") in the 
vicinity of the true value of X+ to get tremen- 
dous saving in the grid points required for a finite 
difference analysis. In so doing the differential 
Helmholtz equation is transformed to a new form, 

Eq. (7), which requires fewer grid points In its 
solution. 

As another example, to show how to pick X+, 
consider the problem of predicting the attenuation 
in a soft wall duct, where rf* - is an unknown. 
Calculations for the attenuation were made in Ref. 

5 using the assumption that p + equals n« Fig- 
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(15) 


ure 5 fro® Ref- 5 shows a comparison of the ana- 
lytical and wave envelope calculation for the 
optimum soft-wall duct attenuation for various 
H and L/H values. 

Excellent agreement between the analytical 
and numerical calculations were obtained for the 
two-dimensional duct example in Ref. 5. For the 
r\ ° 5 and L/H = 6 case in Fig. 5, the conven- 
tional finite difference theory required 3600 grid 
points while the wave envelope difference theory 
required only 100 grid points. Thus, a savings 
of 3500 grid points over the conventional differ- 
ence theory is obtained. At lower n and L/H, 
the savings was smaller. For this type of cal- 
culation for any assumed value of n+, it is nec- 
essary to check for a converged answer by increas- 
ing the number of grid points (m,n) in both x 
and y direction, as was done in Ref. 5. 

Boundary Conditions 

In the transverse direction, the acoustic im- 
pedance at the wall as shown in Fig. 1 is defined 
as the ratio of pressure to the transverse veloc- 
ity. In dimensionless form the impedance ratio 
can be expressed as( 5 ) 

C ’ pc “ " i27rn (3P/3y) (10 


p<2) = o 

Axial Acoustic Power 

The sound power which leaves a duct and 
reaches the far field is related to the axial in- 
tensity. As was shown in Ref. 5, the axial inten- 
sity can be expressed in terms of the transformed 
pressure as 

I = £ [p (1)2 + p (2)2] 



It should be noted that the calculated value of 
intensity will be independent of the assumed value 
of ti + . The total dimensionless acoustic power is 
the integral of the intensity across the test sec- 
tion 

Jq K x *y) d y 07) 

By definition the decrease, in decibels, of 
the acoustic power from 0 to x. can be written as 

AdB = 10 log 10 (f^) (18) 

Difference Equations 


Substituting the expression for the transformed 
pressure, Eq. (5), into Eq. (10) yields 

„ "ISEIE (11) 

3y £ 

which Is the same form as Eq. (10) since the n + 
parameter does not enter into the impedance condi- 
tion. 


The exit condition to be presented now allows 
the numerical solution to be compared to the ana- 
lytical results for wave propagation in an infin- 
itely long duct. The entrance region of length L 
in an infinitely long duct with uniform impedance 
will not have reflections at any position in the 
duct. The wave propagation in the entrance region 
of the infinite duct can be represented with a fin- 
ite length of duct L by choosing the exit imped- 
ance at L so that no reflections occur. Conse- 
quently, if the exit condition in the numerical 
analysis is chosen to eliminate reflections at the 
exit of the duct, the numerical and analytical re- 
sults should be in close agreement. The exit con- 
dition developed in Ref. 5 to meet this condition 
was shown to be 


|| - -mr,p 02) 

Substituting the expression for the transformed 
pressure, Eq. (5), into Eq. (12) yields 

« -127T(n“n + )p (13) 


For the entrance pressure profile, the trans- 
formed pressure is assumed to be 

p(l) - 1 (14) 


The continuous system is now reduced to an 
equivalent lumped-parameter system by means of a 
finite difference approximation. Instead of a 
continuous solution for the pressure, the pressure 
is determined at the isolated grid points shown 
in Fig. 1. The difference equations to follow are 
generated using Crcen*s Theorem. 

Central 


The generalized wave envelope Eq. (7) written 
in finite difference form for the i,j point in 
Fig. 1 is 

© !p i,i-i * [ l + 

- 4 O 2 - ' ^K i 


+ [l - 127rW]p 1+1>j 

, /Ax\2 _ 

+ y p i.j+i ■ ° 


Wall 

The difference equation which applies along 
the upper boundary in Fig. 1 is 

(— Vp + ~fl + i2uAxn + l P, , 

\Ay/ ,m-l 2L J r i-l,m 

r /Ax) 2 (27iAx) 2 (n 2 - H+ 2 ) . 12irnAx 2 l .. 

‘ f + W/ " 2 c Ay 

+ | [l - i2irAxn + ]p i+1>m » 0 (20) 
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(25) 


A similar equation applies along the lower wall. 


{A) [p] - [T] 


Exit 


The difference equation which applies along 
the exit plane in Fig. 1 is 


1 /Ax\ 2 

2 W/ P n,j-1 


+ 


[l + i2irixn + ]p n _ 1>;J - [l +(g) 2 


Corner 


_ (2ttAx) 2(^ - n +2 ) + _ n+) 

+ i2irAxn + P A 
J n i J 


, 1 /M 2 

2 \by) 


n,j+l 


(21) 


The difference equation which applies at the 
upper corner (point n,m in Fig. 1) is given by 

1 (“^Vp . + \ fl + 12TTAxn + lp , 

2 VAy/ *n,m-l 2 L J r n-l,m 

_ l ^ + . ( 2TrAx) 2 (n 2 - n +2 ) 

+ + mix (n . „+j 

c Ay 

+ i27rAxn + ]p ■ 0 (22) 

J 

a similar equation exist for the lower corner. 
Axial Intensity 


In terms of the difference notation, the 
axial intensity as given by Eq. (16) can be ex- 
pressed as 


I 




+ P 




+ 


1 

2irn 



SPd) 

3x 




p(l) 


3p(2) 

3x 



(23) 


The total power across a particular cross section, 
as given by Eq. (17), is written in difference 
notation as 


Ki 


m-l 

+ E 

j=2 


\i + 2 I i,m 


Ay (24) 


By evaluating at the entrance and exit 

positions, taking the'log of their ratio and mul- 
tiplying by 10, as indicated in Eq. (18), the 
sound attenuation for the duct is determined. 


Matrix Solutions 


The collection of the various difference 
equations at each grid point forms a set of simul- 
taneous equations which can be expressed in matrix 
notation as 


where (A) is the known coefficient matrix, CpJ 
is the unknown pressure vector, and [F] is the 
known column vector containing the various initial 
conditions. The matrix is complex. At the present 
time, the most efficient computer routine for 
solving Eq, (25) was written by D. W* Quinn at 
Wright-Patterson Air Force Base, Dayton, Ohio. 
Equation (25) can also be solved by subdividing the 
complex matrix into four real matrices. The real 
matrices can be solved using standard elimination 
techniques. 


Results and Discussion 


Calculations are now presented in this section 
to illustrate the use of the wave envelope concept 
in solving two sample problems in treated duct 
attenuation: the effect of source strength radial 

profile; and the effect of variable axial impedance. 
In the following caleualtions n + will be assumed 
to T). In all calculations, the fact that the 
pressure profile are symmetric about the centerline 
was used to reduce the required number of grid 
points in y-direction. 

Radial Source Profiles 


Because of the complexities of the turbo-fan 
noise source, questions exist at the present time 
on the radial distribution of the acoustic energy 
at the fan source, The attenuation was therefore 
calculated for stepped noise source profiles of the 
form illustrated in Fig. 6. The calculation were 
made for a liner which had been optimized assuming 
a uniform source profile* The sensitivity of the 
liner attenuation for various size of the steps was 
determined at three different frequencies. 

Figure 7(a), (b), and (c) illustrates the re- 
sults for the three cases of n = 1, 2, and 5, 
and all at L/H =1. In each of these figures, the 
abscissa represents the radial position of the step 
in the duct. The attenuation for a duct with a 
uniform radial profile is represented by the far 
left ahnd portion of the curve at y = 0. As the 
radial position of the step moves to the right 
towards the absorbing wall, y > 0, a greater por- 
tion of the energy will be absorbed. As seen in 
Fig. 7(a), for a low frequency source (n = 1) 
a modest (3 dB) increase in attenuation is seen 
for either profile as compared with the attenuation 
of the uniform profile. However, as shown in 
Figs* 7(b) and 7(c), for the smaller wavelength 
sound sources (r\ =» 2 and 5) the attenuation for the 
stepped profiles is much larger than for the uni- 
form profiles. 

Higher order modes have been shown by analysis, 
in general, to have higher damping coefficients 
than the lower order modes. Since such a step 
pressure profile can be considered to be composed 
of various higher-order radial modes, and since 
higher-order modes (stepped towards wall) damp 
more quickly than the lower-order modes, the re- 
sults seen here are in qualitative agreement with 

analytical theories^). 

In a practical problem, if the noise source 
profile of a turbo-fan engine is skewed with a 
higher concentration of acoustic intensity closer 
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to the wall, the attenuation predicted assuming 
a uniform source would under-perdict the measured 
attenuation. 

Wall Impedance Variation 

To eliminate the need for heavy, expensive 
and otherwise undesirable splitter rings and also 
to reduce the required length of wall treatment 
(cowl length), recent research is concerned with 
Increasing the attenuation of vail treatment in a 
fixed length of liner. One practical approach to 
Increasing the attenuation of a liner (?) is to 
subdivide the liner into several different segments 
as shown in Fig. 8. Next, to obtain higher atten- 
uations than are possible with a uniform liner, the 
liner segments are jointly optimized to maximize 
the. noise attenuation. 

In optimizing the stepped liners, the starting 
point was the optimized uniform impedance. The 
first liner segment was then optimized while hold- 
ing the remaining segments at the uniform value. 
Next, holding the first segment at its new value, 
the iteration process was continued on the remain- 
ing segments. The process was repeated for the 
whole liner until the change in the duct attenua- 
tion was less than 5%. 

Figure 9 shows the optimum attenuation deter- 
mined herein for stepped impedances of two and 
three steps for a dimensionless frequency of 
11 *= 1 and a liner length-to-height ratio L/H of 
0.5. As seen in Fig. 9, the two and three stepped 
liners have increased the attenuation by 5 and 6 dB 
respectively. Figure 10 shows the optimum atten- 
uation for a two step duct at the higher dimension 
less frequency of T) ~ 5. These calculations for 
T) = 1 and r\ B 5 show that the stepped impedance 
ducts have theoretical attenuation 2.0 to 2.5 
times the attenuation of a uniform impedance liner 
of the same length. The dotted part of the curve 
shown in Fig. 10 is an numerical instability that 
results from insufficient grid points. As shown 
in Ref. 5, the number of grid points in y-direction 
should be increased for increasing n« Unfortu- 
nately, the present computer program could not 
accomidate the required increase in grid points. 

In both the two and three-stepped liners for 
D «= 1 and a two stepped liner for n - 5, the 
resistive component of impedance optimized near 
zero (0.01) in the first segment of the liner. As 
a result, the first portion of the segmented liner 
removes less acoustic energy than the uniform 
liner. This can be seen in Fig. 9 by observing 
that the curves for both the two and three-step 
liners fall below the curve for the uniform liner 
for x less than 0.3. This is to be expected 
since there is no resistive component of impedance 
to dissipate the energy in the first segment of 
the stepped liners. 

The reactive component of the first section 
of the two and three stepped liners does not differ 
substantially from the optimum reactance for a 
uniform liner with a plane wave input, or from the 
optimum reactance associated with the least atten- 
uated mode in a uniform duct. 

The reasons for the enhanced attenuation for 
the stepped duct liner over the uniform liner can 
be reduced from the work in Ref. 4 concerning the 
acoustic boundary layer and the results presented 


acoustic boundary layer and the results presented 
earlier on stepped noise source profiles. In 
uniform impedance liners, an "acoustic pressure 
diffusion boundary layer" is formed along the walls. 
The acoustic pressure at the centerline of the 
duct is "unaware" of the soft wall until the acous- 
tic diffusion boundary layer at the upper and lower 
walls meet. This is illustrated by the results 
shown in Fig, 11, This phenomenon is analogous to 
the problem of a developing laminar velocity pro- 
file in the entrance region of a pipe. It takes a 
prescribed length down the pipe before the growing 
laminar boundary layers at the walls meet in the 
center of the pipe to establish the fully developed 
parabolic laminar flow profile. Apparently the low 
resistive entrance impedance induces large trans- 
verse movement of the acoustic energy which break 
up the "acoustic diffusion boundary layer" and 
redistributes the acoustic energy towards the wall 
in the entrance section of the duct. The redistri- 
bution of acoustic energy towards the wall is il- 
lustrated for the two-stepped liner in Fig. 12, 
which compares the sound pressure level (SPL) for 
the two-stepped and the uniform impedance liners. 

As seen in Fig. 12, the SPL in the initial segment 
of the duct increases in the vicinity of the wall. 

As was shown in the earlier section of this report 
on radial distribution effects of sound sources, 
if the acoustic pressure profile Is peaked nearer 
to the wall, the acoustic energy in the duct can 
be more effectively absorbed. This phenomenia 
could be called pressure redistribution. The pres- 
sure source is conditioned for better absorption 
by the downstream liner segments. As seen in 
Fig. 12, the peak SPL level which occurred at 
x “ 0.25 is quickly lowered below the uniform SPL 
level at the point x = 0.375. 

To fabricate a low-resistive element, large 
open area in the acoustic liner is required. For 
the no flow case considered here, a closely packed 
array of quarter-wave tubes might be ideal in this 
application since these tubes can be tuned to give 
the proper reactance yet still have a very low re- 
sistance. However, a practical liner may have a 
higher resistance than the desired optimum value. 

In this case, the efficiency of the liner would be 
reduced. The sensitivity of the liner attenuation 
to the entrance resistance is shown in Fig. 13. As 
seen in Fig. 13, the efficiency of the liner de- 
creases with increasing resistance of the first 
segment. Consequently, both theoretical and exper- 
imental work is needed on the design of low-resis- 
tance liners. 

Combined Stepped Segments and Source Profile 

Finally, suppose a stepped noise source pro- 
file of the form shown in Fig, 6 exits in the duct. 
Such a profile might result because more noise is 
generated at the tip of the rotor or because a 
centerbody channels more acoustic energy to the 
vicinity of the wall. Suppose we also want to use 
a two-stepped segmented liner. For a stepped pro- 
file which extends 30£ inward from the wall, the 
optimum impedance was determined for a two stepped 
duct. The results shown in Fig. 14, are compared 
with the uniform impedance optimized for both a 
step input and a plane wave input. In this example, 
the resistance of the initial segment is at 0,11. 
This is higher than the resistance required for the 
plane wave source input for a two-stepped liner. 
Recall that the resistance for the two-step liner 
for a uniform plane wave input was shown to be 
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0.01 In Fig. 9. Most likely, since the acoustic 
pressure is already higher in the vicinity of the 
wall, the optimum resistance in the first section 
of the stepped duct is higher so that acoustic 
energy can be immediately absorbed. Clearly, the 
optimum impedance distribution will be a function 
of the source noise distribution. 

Recall that for a uniform duct with a plane 
wave input, for r\ “ 1 and L/H “0.5 the atten- 
uation is -4 dB, as shown in Fig. 9. Comparing 
this value to the nearly -40 dB attenuation shown 
in Fig* 14 for a stepped noise source profile and 
a segmented liner, the combined consideration of 
both source profiles and segmented liners can be 
a significant in the design of an acoustic liner. 

Conclusions 


Because the wave envelope technique greatly 
reduces the required computer storage and solution 
time, the technique was conveniently applied to 
the acoustic problems of skewed inlet pressure 
waves and ducts with axial variations in impedance. 
For both plane wave uniform noise sources and 
stepped noise sources, segmented liners are shown 
theoretically to increase the attenuation as com- 
pared to a uniform liner by redistributing the 
energy to the wall so that the acoustic energy can 
be more efficiently absorbed by the liner* To 
accomplish this redistribution of energy, the first 
element of the segmented liner must have a value 
of resistance near zero. For higher initial seg- 
ment resistance, the attenuation significantly 
falls below the optimum value. The optimization 
of segmented impedance components and the optimum 
attenuation of the liner are shown to be strong 
functions of the source noise distribution. 
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Dimensionless axial coordinate, x 

Figure Z - Analytical and numerical pressure profiles for one- 
dimensional sound propagation in hard-wall duct. Dimensionless 

frequency y ■ L P ■ + \P^\ 


x + 1 



Figure 3. - Typical pressure profile for sound propagation 
in a two dimensional soft-wall duct for dimensionless 
frequency v = 1 and duct length L/H = 3. 
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Figure 5 . - Effect of axial length and frequency on attenuation at 
optimum impedance in two-dimensional duct assuming 17 + =77 
(ref. 5 ). 
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Figure 6. - Illustration of noise source profile in inlet of turbo-fan engine. 
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t Radial position of step, y Wall 

Duct 

<b> 77 = 2, 0 =0.47, x--l. 32 (plane wave optimum imped- 
ance). 

Figure 7. - Effect of transverse stepped source profiles on 
duct attenuation for L/H = 1 and zero and 0. 3 step for 
m » 14, n = 10 at the f\ and impedance values indi- 
cated in each figure. 
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Axial position, x 

Figure 10. - Impedances optimized for an UH = 0 5 and 
V =5. 
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Figure 11. - Pressure amplitude fluctuations along duct at 
optimum impedance for dimensionless frequency rj = 1, 
dimensionless duct length L/H = 1, number of grid rows 
m = 10, and number of grid columns n » 10 (ref. 5). 
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Figure 12. - Radial SPL profiles for uniform and two-step duct liners 
for ii = 1 and L/H = 0. 5 for m = 5, n = 10. 
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Figure 13. - Effect of initial stepped specific acoustic re- 
sistance 0 ( 1 ) on maximum attenuation in a two-step 
duct for 17 = 1 , L/H « 0. 5, m = 5 (symmetry), n = 10. 
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Figure 14, - Impedances optimized for a stepped noise source profile and two 
segment liner for r\ - 1 and L/H * 0,5. 


